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ABSTRACT 



Context. A significant number of the Wolf-Rayet stars seem to be binary or multiple systems, but the nature of many of them is still 
unknown. Dedicated monitoring of WR stars favours the discovery of new systems. 
Aims. We explore the possibility that WR62a is a binary system. 

Methods. We analysed the spectra of WR62a, obtained between 2002 and 2010, to look for radial- velocity and spectral variations 
that would suggest there is a binary component. We searched for periodicities in the measured radial velocities and determined orbital 
solutions. A period search was also performed on the 'All-Sky Automated Survey" photometry. 

Results. We find that WR62a is a double-lined spectroscopic binary with a WN5 primary star and an 5.5-6 type secondary com- 
ponent in orbit with a period of 9.1447 d. The minimum masses range between 21 and 23 M Q for the WN star and between 39 and 
42 M G for the O-type star, thus indicating that the WN star is less massive than the O-type component. We detect a phase shift in the 
radial-velocity curve of the He n /14686 emission line relative to the other emission line curves. The equivalent width of this emission 
line shows a minimum value when the WN star passes in front of the system. The analysis of the ASAS photometry confirms the 
spectroscopic periodicity, presenting a minimum at the same phase. 

Key words. Stars: binaries: spectroscopic - Stars: individual: WR62a - Stars: Wolf-Rayet 
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Wolf-Rayet stars frequently occur in binary systems. In fact, 
multiplicity seems to be one of the mo st striking properties of 
massive stars (cf. ISana &Evansll20Tlb .l van der Hucht (2001) 
has found 38% binaries and possible b inaries among a total of 
227 WR stars. lLanger & Hegerl (1 19991) argue that most of the 
Supernova Type Ib/c occur in WR + OB interacting binary sys- 
tems, and given the statistics derived from the SN observed, it 
would be expected that more WR stars are detected in binary 
systems than as single stars. Another reason is provided by the 
observation of late-type WC stars (WC8 - 9). These stars exhibit 
thermal emission from hot dust, and it is generally accepted that 
the production of dust is a consequence of the binarity, where 
the dust is created via the collision of hydrogen-ric h and carbon - 
rich winds from the respective OB and WC stars (IUsovlll99l1) . 
Moreover, a systematic high-r esolution spectrosc opic survey of 
O and WN stars (OWN Survey; iBarba etalJ2010h has found that 
nearly 60% of the sample are candidates to belong in binary or 
multiple systems. Therefore, the binary fraction is expected to 
be higher than what was found by van der Hucht, that is, the bi- 
narity appears to play a key role in the formation of WR stars. 
In recent years, the discovery of new binary systems has been 
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favoured thanks to technological advances, allowing us to un- 
derstand a little more about the nature of massive stars. 

The stellar mass is a fundamental astrophysical parameter 
whose knowledge, together with the rate of mass loss, chemical 
composition, and rotation, determines the properties and evolu- 
tion of stars (Mev net & Maederll2005l) . Observing binaries, par- 
ticularly double-lined spectroscopic binary systems (SB2) that 
are eclipsing, is the most direct way to obtain accurate stellar 
masses. 

We are carrying out spectroscopic monitoring of faint south- 
ern galactic Wolf-Rayet stars in order to discover new binary 
systems. One of the most conspicuous radial-velocity (RV) vari- 
ables observed in our sample is WR62a. This star (SMSNPL1 1, 
C2000 = 14 / '32 m 38 s .18; 6mm = -61°29'55'.'4; v=13.8) was dis- 
covered by IShara et ail (1 19991) . They classified it as WN4-5o, 
which means that they did not identify absorption lines in the 
spectrum. 

In this paper, we present the first detailed spectral analysis 
for WR62a. We demonstrate that it is a short-period SB2. We 
also report a phase shift in the orbital solution of one of its emis- 
sion lines and discuss the photometric variability of the star. 



2. Observations and data reduction 

We have obtained a total of 28 spectral images of WR62a. The 
images were acquired with the 4-m V. Blanco telescope at Cerro 
Tololo Inter- American Observatory (CTIO), Chile; the 2.5-m du 
Pont telescope at Las Campanas Observatory (LCO), Chile; and 
the 2.15-m J. Sahade telescope at Complejo Astronomico El 
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Table 1. Details of the observations. 



Date-Obs. 


n 


Sp. coverage 


Observat. 


Telesc. 


Spectr." 


Grating 


Detector 


Dispersion 


Resolution* 


UT 




[A] 




[m] 




[lmirr 1 ] 




[A pix" 1 ] 




2002, Apr. 22 


1 


3940-5600 


CASLEO 


2.15 


REOSC 


600 


Tekl024 


1.63 


1000 


2007, Mar. 29 - Apr. 3 


6 


3650-6700 


CTIO 


4 


R-C 


632 


Loral 3k 


1.01 


1300 


2008, Apr. 19-23 


6 


3650-6700 


CTIO 


4 


R-C 


632 


Loral 3k 


1.01 


1300 


2009, Mar. 24-27 


4 


4030-5590 


CASLEO 


2.15 


REOSC 


600 


Tekl024 


1.63 


1000 


2009, Apr. 20-21 


2 


3800-5530 


LCO 


2.5 


B-C 


1200 


Marconi#l 


0.79 


2500 


2009, Jul. 19-27 


4 


3520-5250 


LCO 


2.5 


B-C 


1200 


Marconi* 1 


0.79 


2500 


2009, Aug. 14-16 


3 


4040-5700 


CASLEO 


2.15 


REOSC 


600 


Tekl024 


1.63 


1000 


2010, Apr. 10-13 


2 


3930-5600 


CASLEO 


2.15 


REOSC 


600 


Tekl024 


1.63 


1000 



n : Number of spectra obtained per run. 

a : Details of the spectrographs can be found in the User's Manuals of the respective Observatories. 

b : The spectral resolutions (R=A/AA) were measured using AA as the FWHM of the calibration lamp emission lines. 



Leoncito (CASLEO fl Argentina, between 2002 and 2010. See 
Table Q] for a summary of the observations. We have acquired 
spectra at CTIO using the Ritchey-Chretien (R-C) spectrograph, 
at CASLEO employing the REOSC spectrograph, and at LCO 
using a Boiler & Chivens (B-C) spectrograph. The pro cedure for 
obtain ing LCO spectra are also described in detail by ISota et aD 
(1201 lh . Slight differences in the spectral ranges listed in Table Q] 
are due to changes in the grating angle among the observing 
runs. 

Typical exposure times for spectra were between 15 and 40 
min, resulting in spectra with signal-to-noise ratios S/N ~ 50 
- 150. Comparison lamp spectra of He-Ne-Ar (or Cu-Ar with 
REOSC spectrograph) were observed immediately after or be- 
fore, at the same telescope position as the stellar exposures. Bias 
and flat-field images were also obtained each night. We did not 
observe a telluric standard star since there are not many atmo- 
spheric absorption lines in the wavelength range of our spectra. 
All spectra were processed with IRAF0routines. 



3. Results and discussion 

3.1. The spectrum of WR 62a 

Our spectra reveal that WR62a is a double-lined WR+OB bi- 
nary system. Several absorption lines are detected in our spectra. 
The spectra obtained at different epochs show that the absorption 
lines have different radial-velocity shifts, from those of the emis- 
sion lines, suggesting that the system is an SB2 . We used an iter- 
ative p rocedure similar to the one developed bv lMarchenko et al] 
( 1998) to separate the composite spectrum into the spectra from 
the two components of the binary. First, we begin by shifting 
the spectra to its respective radial velocity and calculating the 
mean spectrum to be used as a template of the WN. Second, 
the RV-shifted mean WR spectrum is subtracted from the origi- 
nal spectra, obtaining the remaining O-type spectra. Then, they 
are shifted to their corresponding RV of the absorption lines and 
combined to get the template of the O-type star. After that, the 
(RV-shifted) mean O-type spectrum is subtracted from the ob- 
served spectra and a new purer WN-type spectrum is obtained. 
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That is to say, the spectra are shifted, combined, and subtracted 
appropriately to compute the individual spectra of the stellar 
components. We noted that the method converges after a few 
iterations. 

After separating the respective spectra (see Fig. [TJ, we per- 
formed a more detailed spectral classification. Relative inten- 
sities of N m, N iv, He n, and C iv emission lines in the WR 
spectrum indic ate a WN5 spect ral type (according to the cri- 
teria given by Smith et al. 1996), confirming the classification 
proposed bv lShara et al. 1 (119991) . 

In the disentangled spectrum of the O-type star, we identi- 
fied absorption lines of Ho-, H/3, Hy, Hc5, He i AA 4026, 447 1 , and 
5875, and also Hen AA 4200, 4542, 4686, and 5411, as well as 
the C in 5696 emission line. The comparison between the He i 
/14471 and the Hen A4542 absorption lines indicates that the 
O-type star is earlier than 07 spectral type; i.e., the intensity 
of the former He i line is fainter than that of the He n line. An 
improved spectral classification was achieved through the com- 
parison of the O star spectrum with the MK standard s from the 
new Atlas for spectral classification (ISota et al.ll2.01 ll) . Thus, an 
O 5.5-6 type is determined. A luminosity class is not possible to 
determine owing to the profile variations in the emission lines of 
the WN-type component which introduce residuals in the result- 
ing O-type spectrum (for example around 4600 - 4700A). 



3.2. Analysis of radial velocities 

We have determined the RVs of the WN component of WR62a 
by measuring the position of the stronger emission lines, i.e. N iv 
,14058, N v ,14604, and He n ,1,14200, 4542, 4686, and 541 1 . For 
the O-type star, we determined a mean RV value for each spec- 
trum by averaging the RV measurements derived from the ab- 
sorption lines. These measurements are listed in Table [2] 

The strongest emission line in our spectra corresponds to 
He ii ,14686, so we could measure this line even in spectra with 
poor signal-to-noise ratios. The RVs of this line show large vari- 
ations from night to night, suggesting an orbital period of a few 
days. We have searc hed for periodicities by me ans of the al- 
gorithm published by iMarrac o & Muzzio (1980). This method 
divides the observations into several phase intervals and com- 
putes the variance with respect to the best-fitting straight line in 
each interval. The most probable period obtained corresponds to 
9.144 d (Fig. [2]>, and it is used as ian initial value to feed orbital 
solutions for the RV data. 
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Fig. 1. (a) Normalized spectrum of WR62a before the "disentangling" method, (b) the WR star component after removing the 
O-type star, (c) the remaining spectrum of the O-type star. The spectra are shifted in the intensity axis for a better comparison. 




Period [d] 

Fig. 2. Pe riodogram obtained for the He n /14686 emission line 
with the M arraco & Muzziol (11980b code. The dispersion (Y- 
axis) is the parameter used to estimate the quality of the fitting 
for each trial period. This value is the total variance of the best- 
fitting straight line to the observations in each phase interval. 



We performed orbital fits to each dataset of RVs, using the 
GBART0 program. It is clear that the determined RVs for the differ- 
ent spectral features show the same periodicity (see Table |3), and 



thus, we adopted the straight mean of the periods as the definitive 
value for the orbital period. The ephemeris adopted for WR62a 
is 



T = 2,455,038.9 + 9.1447dx£ 



(1) 



where Tq is the time when the WN star passes in front of the sys- 
tem, as in the normal convention for eclipsing WR+OB systems. 
Then, we calculated the orbital solutions with the period fixed to 
this average. The orbital parameters derived for each RV dataset 
are presented in Table [4] (with their corresponding probable er- 
rors). To determine the To of the system, we neglected the very 
low eccentricity, averaging all the TRVmax values of the emission 
lines of Tableland calculating the Tq as < TRVmax > +0.75 x P, 
which coincides with the primary minimum of the lightcurve (as 
will be shown in Sect. l3.41 >. Some orbital solutions are illustrated 
in Fig. |3] The RV curves are phased with the ephemeris from 
Eq. 1. 

In Fig.[3]and Table|4] a phase shift between the Hen /14686 
line and other emission lines is noticeable. As a result, the 
T^RVmax value for this line was excluded from the calculus of the 
ephemeris. This phase-shift effect is also det ected in other close 
WR binary systems e.g. HD 90657 (WR 21: INiemela & M offat 
11982k iGamenl 12004 . HP 94546 (WR 31; lGamenlj2004 . and 

04), 



ii^oza ivjui iicn z,uutj, i iu y^tj^tu yvyjv ji, vj allien z.uu^, anu 

HDE 320102 (WR 97; INiemela et al]fl995l: iGamerEool . and 



3 gbart is an improved version of the program originally written by 
iBertiau & Grobbenl ([1968) and developed by Federico Bareilles. 



usually interpreted as a sign of colliding winds. The amplitudes 
and systemic velocities of the orbital motion determined from 
the RV of different emission lines show differences. The ampli- 
tude of N v i4944 shows the lowest value, while Hen /Lt4200, 
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Table 2. Observed heliocentric radial velocities of WR62a. The mean FWHM and mean EW of each emission line measured from 
the separated spectra are indicated in the last two rows, respectively. 



T-TTFl 
tl.1 1 J 


Phase 


"NT T17 

IN IV 


IN V 


NT \- 
IN V 


V, IV 


MC II 


rlc II 


rlc II 


IT, , Try 

rlc II 


Mean Abs [nj 


2,450,000+ 




4057.66 A 


4603.73 A 


Af\ A A C f A 

4944.56 A 


C O 1 1 Ok O A 

5811.98 A 


A f O C /" O A 

4685.68 A 


5411.52 A 


A 1 t~\f~\ O ^7 A 

4199.87 A 


A C A 1 CC\ A 

4541.59 A 








n —In 

[km s 1 ] 


rl - 1 l 

[km s '] 


rl -It 

[km s '] 


n —In 

[km s 1 ] 


rl —In 

[km s 1 ] 


rl -It 

[km s 1 ] 


n —It 

[km s '] 


n —It 

[km s '] 


n -It 

[km s '] 


2386.72 


0.98 


-213 


49 






113 


55 






-43 ± 27 (5) 


4188.89 


0.05 


-65 




-41 


59 


254 


235 






-113 ±22 (8) 


4189.82 


0.15 


53 


220 


106 


155 


278 


306 


308 


387 


-204 ± 32 (10) 


4191.84 


0.37 


-8 


213 


96 


185 


134 


223 


297 


397 


-175 ±29 (10) 


4192.85 


0.48 


-161 


9 


-8 


-16 


-61 


86 


162 


261 


-83 ± 22 (9) 


4193.65 


0.57 


-275 


-119 


-199 


-155 


-162 


-62 


42 


-30 


19 ± 59 (9) 


4193.87 


0.59 


-228 


-143 


-280 


-135 


-250 


-192 


-72 




18 ± 52 (9) 


4575.84 


0.36 


87 


187 


102 


182 


133 


259 


374 


359 


-157 ±36 (10) 


4576.82 


0.47 


-122 


48 


-32 


-19 


-53 


79 


90 




-85 ± 28 (6) 


4577.73 


0.57 


-285 


-116 


-193 


-178 


-199 


-119 




-14 


19 ± 38 (9) 


4578.80 


0.69 


-403 


-254 


-346 


-257 


-282 


-284 








4579.59 


0.77 


-379 


-255 


-276 


-257 


-191 


-234 


-246 


-215 


31 ±37 (9) 


4579.85 


0.80 




-257 




-290 


-173 


-250 


-262 


-230 


36 ± 35 (9) 


4914.83 


0.43 


-8 


173 


101 




49 


130 




307 


-75 ± 26 (5) 


4915.82 


0.54 


-134 


-14 


-161 




-89 


-25 


114 


139 


25 ± 18 (3) 


4916.78 


0.65 


-255 








-251 


-208 


-181 


-86 


98 ± 14 (4) 


4917.74 


0.75 


-406 


-216 


-343 




-183 


-282 




-212 


91 ± 23 (5) 


4941.57 


0.36 


91 


254 


150 




218 


278 


375 


362 


-151 ±24(7) 


4942.71 


0.48 


-163 


19 


-97 




-12 


-21 


166 




-77 ± 28 (8) 


JUjl .04 


U.Z I 


SI 1 
O I 


Z / J) 


1 J J 




3k) 1 


3 no 

JUZ 


JiZO 




— ZU4 ± ZI \Q) 


5036.59 


0.75 


-372 




-288 


-212 


-168 


-227 


-249 


-205 


102 ± 24 (7) 


5038.58 


0.97 


-204 


-41 


-167 


-185 


131 


33 


-32 


109 


-97 ± 36 (3) 


5039.61 


0.08 


-31 


123 


-1 




263 


218 


228 


244 


-168 ±33 (7) 


5057.52 


0.04 




170 


-118 




352 


314 






-137 ±40 (7) 


5059.52 


0.25 




287 


175 




356 








-183 ± 28 (6) 


5060.49 


0.36 




191 


60 




187 


237 






-151 ± 17 (7) 


5296.79 


0.20 


153 








267 






373 


-195 ± 33 (6) 


5299.66 


0.51 










-113 








31 ± 1 (4) 


FWHM [A] 




18.5 


10.4 


19.7 


33.9 


24.7 


22.4 


17.3 


18.7 




EW [A] 




-4.6 


-1.6 


-1.5 


-2.2 


-26.1 


-3.8 


-1.2 


-1.6 





a: the errors are the standard deviation of the mean, and the (n) values indicate the number of absorption lines measured. 



Table 3. Periodicities derived from the radial velocity of N iv, 
N v, C iv, He ii emission lines and the mean absorption. 



Feature 


P[d] 


error [d] 


Niv 44058 


9.1444 


0.0005 


Nv 44604 


9.1455 


0.0005 


Nv 44944 


9.1443 


0.0009 


Civ 45812 


9.146 


0.002 


Hen 44686 


9.1447 


0.0005 


Hen 45411 


9.1447 


0.0003 


He ii 44200 


9.145 


0.001 


Hen 44542 


9.140 


0.002 


Mean Abs. 


9.1415 


0.0005 



4541, and 541 1 emission lines appear somewhat larger. The sys- 
temic velocities of emission lines (except N iv /14058 and N v 
/14944) are redshifted relative to that of the absorption lines. This 
is a known phenomenon in WR+OB systems, which indicates 
that emission l ines are being formed in different regions o f the 
WN wind (cf. iNiemela & Moffad[l982t INiemela et alj[l995h . or 
perhaps are still affected by the absorption lines of the secondary. 

From Tableland Fig.[4]it is also evident that the absorption 
lines (the mean of H, Hei, Hen) move in anti-phase with the 
emissions, thus they belong to a secondary component. Besides, 
the semi-amplitude of the orbit of this companion is lower than 
the WN one (regardless of emission line used), which means that 
the secondary is more massive than the WN component. 

We performed fits of the orbital solution for both compo- 
nents together, by means the gbart program. It is widely known 



that the RVs of the different emission lines do not strictly follow 
the same orbital motion in WR binary systems, which is often 
attributed to their different line-forming regions. We assumed 
that the orbital motion of the WN star is better represented by 
the RVs of the highest ionization emission lines, namely, N iv A 
4058 A and N v A 4604 A. The O-type star is depicted by the RVs 
of the mean of some absorption lines. To perform the simultane- 
ous fit of both RV curves, we shifted the emission-line RVs to 
match the systemic velocity derived for the O-type component. 
Therefore, the RVs of Niv and Nv were shifted by 80 kms" 1 
and -86 kms -1 , respectively (see Table |4j. These solutions are 
shown in Table [5] and Fig. [4] Thus, the secondary component is 
confirmed as the more massive star. Depending on which orbital 
solution is considered, the mass of the O-type star ranges be- 
tween 39 and 42 M , and the WN mass, between 21 and 23 M & . 
The mass ratio (q) obtained agrees with the valu es determined 
in other WN5+OB systems (Ivan der Huchd 120011) . The low ec- 
centricity determined using different datasets indicates that the 
orbit is almost circular, although such a low value should not be 
neglected. 

3.3. Analysis of spectral features 

In Fig. [5] we show the Hen /14686, line at four different orbital 
phases; the absorption line of the secondary can be seen to be 
moving in anti-phase with respect to the emission line peak. In 
some phases, this emission line shows an asymmetric profile, 
perhaps indicating some additional contribution from the wind- 
wind colliding region. We measured the equivalent width (EW) 



Collado et al.: The new Wolf-Rayet binary system WR62a 5 



Table 4. The orbital parameters (with their respective probable errors) derived from the radial velocity curves of N iv, N v, C iv, 
He ii emission lines, and the mean absorption. 



Parameter 


Niv 


Nv 


Nv 


Civ 


Hen 


He ii 


Hen 


Hen 


Mean Abs. 




4058 


4604 


4944 


5812 


4686 


5411 


4200 


4542 




P[d] 










9. 1447 ±0.0014 










V [fans- 1 ] 


-145 ±5 


21 ±5 


-94 ±5 


-36 ±7 


42 ±4 


47 ±4 


80 ±6 


133 ±6 


-65 ±3 


Kfkms- 1 ] 


261 ±7 


267 ±7 


246 ±7 


275 ± 13 


272 ±6 


303 ± 5 


321 ± 8 


320 ±8 


143 ±5 


e 


0.07 ± 0.02 


0.05 ± 0.03 


0.08 ± 0.02 


0.11 ±0.04 


0.07 ± 0.02 


0.09 ± 0.02 


0.08 ± 0.02 


0.11 ±0.03 


0.08 ± 0.04 


a) [deg] 


12 ±22 


160 ± 29 


24 ±20 


-14 ± 17 


159 ± 16 


183 ± 11 


166 ± 19 


160 ± 14 


319 ±21 


TRVmax MP 


5032.1 ±0.5 


5032.0 ± 0.7 


5032.2 ± 0.5 


5032.1 ±0.4 


5031.4 ±0.4 


5031.8 ±0.3 


5032.3 ± 0.5 


5032.2 ± 0.4 


5036.1 ±0.5 



* TRVmax = Time of maximum RV in Heliocentric Julian days -2, 450, 000 
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Fig. 3. Radial velocities of N iv ,14058 , N v ,14604, C iv A5S12, 
and Hen /14686 emission lines phased with the ephemeris of 
Eq. 1. Continuous curves represent the orbital solutions from 
Table|4] The phase shift of the Hen line is noticeable. 



and full-width at half-maximum (FWHM) of the Hen A4686, 
N iv ,14058, and N v ,14944 emission lines. We detected a certain 
variability of these line parameters, which seems to be modu- 
lated through the orbit (see Fig. |6). The behaviour of the EW of 
the He n line is very noticeable: it decreases around phase 0=0.0, 
when the WN star is in front of the system. This phenomenon 
is also found in other WR+OB binary systems, i.e. WR2 1 and 
WR47 (lGamenll2004l) . and HD 5980 dMoffat et alJll998h . ex- 
plained as if the emission line had an additional non-stellar 
component, originating in the colliding-winds region, which is 
eclipsed during the conjunction phase. 
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Fig. 4. The RV curves for both components in WR62a. The open 
circles show the RVs for N iv A4058 A and N v A 4604 A emis- 
sion lines. The filled circles show the mean RVs of some absorp- 
tion lines for the secondary one. RVs of the emission lines were 
shifted to match the systemic velocity derived from the absorp- 
tion lines alone (see text). The ephemeris given in Eq.l was used 
to generate the RV curves. 

Table 5. Orbital solutions for both component of WR 62a 





Niv+abs. 


N v+abs. 


P[d] 


9.1447 (fixed) 


e 


0.05 ± 0.02 


0.05 ± 0.02 


K WR [fan s~'] 


257 ±7 


266 ±6 


K [fans" 1 ] 


143 ±6 


143 ±6 


V [fans- 1 ] 


-62 ±3 


-62 ±3 


m [deg] 


39 ±25 


149 ± 20 


TVeriasI M 


5, 032.9 ±0.6 


5, 035.7 ±0.5 


TRVmax M* 


5, 032.0 ±0.6 


5,031.9 ±0.5 


iWR sin j [R a ] 


46.1 ± 1.1 


47.8 ± 1.1 


a sin i [R Q ] 


25.6 ± 1.1 


25.7 ± 1.1 


M W r sin 3 i [M ] 


21.5±4.8 


22.6 ±5.0 


Mo sin 3 (' [M Q ] 


38.7 ±5.2 


42.0 ±5.0 


q(WR/0 


0.56 ± 0.04 


0.54 ± 0.04 


* HJD-2,450,000 



3.4. Analysis of the available photometry 

As the SB2 orbital solution for WR62a suggests high minimum 
masses for both components of the system and the EW of the 
Hen /14686 emission line has a minimum at one of the con- 
junctions, we would expect to observe photometric variability 
(eclipses, for instance). We, therefore, analysed the data pub- 
lished in the "All-Sky Automated Survey" (AS AS; Poimariski 
1200 lh . which gives the V magnitudes for this star. Again, using 
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Fig. 5. Continuum rectified spectra of the He n /14686 emission 
line with the superimposed absorption observed during different 
orbital phases of the binary. A noticeable anti-phased movement 
of absorption and emission features is shown. 
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Fig. 6. EW and FWHM measured in the Hen ,14686 (filled cir- 
cles), Niv /14058 (crosses), and Nv ,14944 (open circles) emis- 
sion lines. 



the code of iMarraco & Muzziol ([T980), we searched the AS AS 
data for periodicities and obtained the most probable period of 
9.1438+0.0009 d, very similar to the period obtained from our 
RV analysis (see the periodogram in Fig. |7). In Fig. [8] we have 
plotted the AS AS V magnitudes against the orbital phases gener- 
ated using the ephemeris given in Eq. 1 . We see a V-shaped dip at 
<p = 0.0. A much shallower dip is marginally detected at <p — 0.5. 
To improve the visualization of the variations, we also show the 
data points averaged in phase bins of 0.01 in Fig. [8] The shape of 
the dip at - 0.0 calls the non-phot ospheric eclipsing binaries 
analysed by Lamontagne et al.l dl996l) . In that paper, the authors 
analysed the eclipsing binary systems WR21 and WR47, where, 
as we described in the previous sections, the He n ,14686 emis- 
sion line also shows a phase shift in RVs and a minimum of the 
EW when the WN star passes in the front of the system. The sce- 
nario where an additional emission component to the He n ,14686 
emission is formed in a colliding-wind region, which is being 
eclipsed by the WN component, should also apply. However, the 
accuracy of the ASAS data does not allow us to determine the 
inclination of the system, hence the absolute masses for both 
components. The V magnitude of WR62a is close to the faint 
limit of the ASAS survey. More accurate photometry is required 
to confirm or rule out the eclipsing nature of this binary. 
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Period [d] 



Fig. 7. Periodogram obtai ned from the ASAS data with the 
Marr aco & Muzziol d!980h code. The most probable period 
found is 9.1438+0.0009 d. 




Orbital phase 

Fig. 8. The lightcurve of WR62a using the ASAS data. The 
individual data-points are shown by crosses; the filled circles 
represent the data-points averaged in phase bins of 0.01. The 
ephemeris derived from our spectroscopic data (Eq. 1) is used to 
calculate the orbital phases. 



4. Conclusions 

From our spectroscopic observations and analysis of the RV of 
the lines detected in the spectra of WR62a, we discovered that 
it is a double-lined binary system with a 9.1447 d orbital pe- 
riod. The relative intensities of the emission lines detected in our 
spectra confirm that the spectral classification of the WN star as 
WN5, while the absorption lines identified for the secondary star 
corresponds to an O 5.5-6 spectral-type. The emission lines of N, 
C, and He ions originated in the WR star follow the same move- 
ment and orbital period, but we detect a slight phase shift in the 
He ii ,14686 emission line, which may be due to asymmetries in 
the line-forming region. 

We performed fits of the orbital solution for both components 
by considering two possible RV datasets of the WN star. From 
these solutions, we confirm that the O-type star is more massive 
than the WN star. The minimum mass derived for the O-type star 
ranges from 39 to 42 M Q , and for the WN-type star from 21 to 
23 M Q , depending on which emission line is used to represent 
the orbital motion of the WN star, i.e. N iv ,14058 orNv/1 4604 
respectively. 
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We detected variations in the EW of the He n A4686 emission 
line with the orbital period, showing a ~30% "dip" around (f> = 0, 
where the WN star is in front. This effect suggests that an addi- 
tional emission component to this line, originated in a colliding- 
wind zone, is being eclipsed. Analysing photometric data from 
ASAS, we also identified a minimum when the WN star is in 
front of the system. Additionally, a very marginal secondary dip 
is detected in the opposite orbital phase (around (f> = 0.5). The 
quality of ASAS data does not allow us to determine the incli- 
nation of the system hence the absolute masses for both compo- 
nents. To confirm the eclipsing nature of WR62a, we need more 
accurate photometry. 

The spectrum of WR62a shows a remarkable resemblance 
to the WR+O binary system WR21 dGamenll2004l) . Similar phe- 
nomena are found in both binary systems. These systems could 
be used as test benches for geometric colliding-wind models. 
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